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The competition between spin-triplet and singlet pairings is studied theoretically for the tight- 
binding Q-/3 bands in Sr2Ru04, which arise from two sets of quasi-one dimensional Fermi surfaces. 
Using multiband FLEX approximation, where we incorporate an anisotropy in the spin fluctuations 
as suggested from experiments, we show that (i) the triplet can dominate over the singlet (which 
turns out to be extended s), and (ii) the triplet gap function optimized in the Eliashberg equation 
has an unusual, very non-sinusoidal form, whose time-reversal-broken combination exhibits a crib- 
shaped amplitude with dips. 



PACS numbers: 74.20-z, 74.20Mn 



Spin-triplet superconductivity is of great conceptual 
importance. It is interesting to ask how the spins in a 
Cooper pair can align, since usually the singlet pairing 
is favored so that some special mechanism should be en- 
visaged to account for a triplet superconductivity. The 
p-wave pairing in supcrfluid 3 He is an outstanding ex- 
ample, where a clear picture of the hard-core interaction 
favoring the triplet exists. 

In the past several years, Sr 2 RuC>J has attracted much 
attention as a strong candidate for triplet superconduc- 
tivity. In a seminal paper, Rice and Sigrist suggested a 
mechanism for triplet pairing in .this material, which they 
call an 'electronic version of He'.H In their scenario the or- 
bital degeneracy causes ferromagnetic spin fluctuations, 
which is considered to favor the triplet pairing. Suh= 
sequent experiments indeed suggested triplet pairingJ3'cl 
However, a new puzzle arose when the spin fluctuation in 
Sr2RuC>4 was found to be antiferromagnetic rather than 
ferromagnetic in a neutron scattering experiments Usual 
wisdom dictates that antifcrrcojiagnetic spin fluctuations 
lead to singlet d-wave pairing.B 

Recently, Kuwabara and one of the present authors,Q 
and independently Sato and Kohmotoja have proposed 
that anisotropy in the spin-pfluctuation, observed in NMR 
experiments for Sr2Ru04,u may lead to triplet p-wave 
pairing. However, simple functional forms for triplet 
and singlet gap—functions were assumed in ref. ^, i.e., 
sinfe a: (-|-jsinfcy)tj for p-wave and cosfc^ — cosk y for 
d x 2 _ y 2 -wave. In ref. ^ the form of the gap when the 
quasi- ID Fermi surfaces are hybridized was discussed 
only qualitatively. The form of the gap function is cru- 
cial in discussing the triplet-singlet competition, so the 
functional form should be optimized. The form of the 
gap function is also important in comparing with exper- 
imental results because, e.g., some recent experiments 
have suggested presence of nodes in the superconduct- 
ing gapj!3l!£l while a thermal conductivity measurement 
suggests an isotropic gap.t 3 ] 

The purpose of the present paper is to determine 
the functional form of the superconducting gap func- 



tion self-consistently for the spin-anisotropy mecha- 
nism mentioned above. We adopt the fluctuation ex- 
change (FLEX) approximationli-3, in which we incorpo- 
rate anisotropy in the spin fluctuations. The FLEX re- 
sults are then fed into the linearized Eliashberg equation 
to obtain the gap function. The optimized triplet gap 
function turns out to have unexpected, non-sinusoidal 
forms, which are in sharp contrast with the k x or sin k x 
gap functions assumed previously. As a result, the am- 
plitude of the gap has a shape of a crib along the Fermi 
surface, which may resolve the controversial experiments. 
The origin of the peculiar behavior of the gap function 
is traced back to the singular fc-dependence in the spin 
susceptibility due to a nesting of two sets of quasi-one 
dimensional Fermi surfaces. 

The ruthenate is essentially a three-band system, 
which arises from d xz orbitals aligned linearly along the 
x axis, d yz along the y, and d xy in the xy plane. The for- 
mer two give rise to the a, (3 bands, which are well nested 
due to the quasi-one dimensionality and causes antiferro- 
magnetic spin fluctuations. In this paper, we concentrate 
on the a, (3 bands, namely, we consider a tight-binding 
model, 
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on a square lattice. Here c™ creates an electron at 
d m (m = xz or yz) orbital, the nearest-neighbor hop- 
ping integral t is along the x{y) direction for d xz (d yz ) 
orbitals. We take t — 1 as a unit of energy. 

We have also included the next nearest-neighbor hop- 
ping t' which corresponds to a weak hybridization. When 
t' 5^ the two sets of quasi-one-dimensional bands anti- 
cross, and two two-dimensional (rounded-square) bands 
result, which are the a and (3 bands. The on-site repulsive 
interaction, U, is considered within each orbitals, and in- 
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terorbital interactions are neglected for simplicityJli The 
band filling is n = 4/3 electrons per orbital in S^RuGi. 

In treating the interaction, we employ the FLEX ap- 
proximation. This method is a kind of self-consistent 
random-phase approximation (RPA) where the dressed 
Green's function is used in.ihp-.RPA diagrams. In the 
multiband version of FLEXjiMZl the Green' function G, 
the susceptibility \, the self-energy £, and the supercon- 
ducting gap function all become 2x2 matrices, e.g., 
Gj m (k, ie n ), where l,m denote d yz or d xz orbitals. The 
orbital-indexed matrices for Green's function and the gap 
functions can be converted into band-indexed ones with 
a unitary transformation. As for the spin susceptibility, 
we diagonalize the 2x2 matrix xtp an d concentrate on 
the larger eigenvalue, denoted as x zz - 

The actual calculation proceeds as follows: 

(i) Dyson's equation is solved to obtain the renormalized 
Green's function G(k), where k is a shorthand for the 
wave vector k and the Matsubara frequency, ie n , 

(ii) -|The fluctuation-exchange interaction V^(q) is given 
asM 
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The effective interactions due to longitudinal (zz) and 
transverse (H — ) spin fluctuations (sp) and that due to 
charge fluctuations (ch) have the forms V s z p z = U 2 \l*, 
V+~ = U 2 xtf-, and V ch = U 2 Xch, respectively, where 
the spin and the charge susceptibilities are 



x z s ;(q) = x i "m-Ux irr (q)}- 1 , 

xt(q)^a X i "(q)[l-Uax i "(q)}-\ 
xM = X i "(q)[l + U X i "(q)}- 1 , 
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in terms of the irreducible susceptibility x irr (g) = 
— -^^ fe G(fc + q)G(k) (iV:number of fc-point meshes). 
Here we have taken account of the anisotropy in the spin 
fluctuation by introducing a phenomelogical parameter 
a a la Kuwabara-Ogata.Q From the NMR experiments, El 
we assume a < 1. 

(iii) then brings about the self-energy, S(fe) = 

jj ~^2 q G(k—q)V^ 1 '(q), which is fed back to Dyson's equa- 
tion, and the self-consistent iterations are repeated until 
convergence is attained. We take 64 x 64 /c-point meshes 
and up to 16384 Matsubara frequencies in order to ensure 
convergence at low temperatures. 

We determine T c as the temperature at which the 
eigenvalue A of the Eliashberg equation, 
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for singlet pairing, 
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for triplet pairing with S z = ±1 (d A- z), and 
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for triplet pairing with S z = (d \\ z). Here d is the 
d-vector characterizing the triplet pairing gap function. 

When the electron-electron repulsion, which causes 
fluctuations, is short-ranged (as for the Hubbard U in- 
teraction) the spin fluctuations are much-stronger than 
the charge fluctuations, i.e., (V sp 3> Vch).cJ When fluctu- 
ations for a certain ('nesting') wave vector Q are pro- 
nounced, the main contribution in the summation in 
eq.(||) comes from those satisfying k — k' ~ Q, which 
should be the case when the Fermi surface is nested. The 
present Fermi surface is indeed well nested due to the 
quasi-one-dimensionality. 

Now we turn to the results summarized in Fig.pl for 
t' = 0.3 and U = 5. In Fig.|(a), the ridges in \G\* de- 
lineate a and (3 Fermi surfaces. These quasi-lD surfaces 
are strongly nested at q = (2ir/3,q y ) and q = (q x , 2n/3) 
(mod 2tt), so that the spin susceptibility is ridged in a 
crib shape as shown in Fig.[j](b), with peaks at the cor- 
ners, q = (±27r/3, ±27r/3). [This is consistent with neu- 
tron scattering experiments. El 

The triplet and singlet gap functions obtained by solv- 
ing the Eliashberg equation are shown in Fig.^(c) and (d) , 
respectively. Remarkably, the triplet gap function takes a 
strange shape: although the symmetry is consistent with 
p-wave, its form is far from sin k x along the rounded- 
square Fermi surface. Rather, it has an almost constant 
amplitude on a pair of parallel sides, k x = ±27r/3 (mod 
27r), of the square with a vanishing amplitude on the 
other pair (k y = ±27r/3) of parallel sides. This applies 
to each of the a and f3 bands. Of course the symmetry 
dictates that the other solution (p y ), rotated by 90 de- 
grees from what is described here (p x ), enters on an equal 
footing as we shall discuss below. 

Why do we have this peculiar behavior for the 
FLEX+Eliashberg optimized gap function ? To begin 
with, superconductivity arises due to pair scattering from 
(k, — k) to (k', — k') mediated by the pairing interaction 
y( 2 )(q), where q = k k' is the momentum transfer. 
From the BCS gap equation we can see that supercon- 
ductivity arises if the quantity 
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reaches unity. Here the pairing interaction 
by 



(2) 



is given is positive and large, where we denote the gap func- 
tion as 4>^ (/i — s for singlet and t for triplet pairing). 



As discussed by Kuwabara and OgataQ, and indepen- 
dently by Sato and Kohmotou, the p-wave pairing, with 
</>t(k)0 t (k + Q) < 0, is favored when spin anisotropy is so 
strong as to realize F S * Z (Q) > 2V+~(Q), i.e., V^Q) > 0. 
However, the present self-consistent calculation shows 
that the situation is a little more involved. A key fac- 
tor is the spin fluctuation that is enhanced along a line 
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FIG. 1. The contour plot of the FLEX result for Green's 
function (|G| 2 ; a), spin susceptibility (b), optimized gap func- 
tion for the spin triplet (<j>t; c) or singlet (4> a ; d) pairing for 
U = 5, t' = 0.3, a = 0.8, and T = 0.02. The left(right) 
panel for the a(/3) band, and white(black) corresponds to pos- 
itive (negative) amplitude in (c,d). 



q = (q x ,2n/3). For a given nodal line (vertical for p x 
pairing), the pair scatterings across one pair of parallel 
sides of the Fermi surface become all favorable as in- 
dicated by O m Fig.||(a). This gives rise to the near- 
constant gap function on that pair of sides of the Fermi 
surface. By contrast, the pair scatterings across the other 
pair of parallel sides (x in Fig.|](a)) lead to 
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in eq.(|lC|) when k x ,k' x have the same sign. Since —V t K 
is negative, the p^-wave gap functions on these sides of 
the Fermi surface is unfavored. This explains the gap 
function shown in Fig.[l](c). 

Now, let us move on to the competing superconduct- 
ing state in the singlet channel. Although Kuwabara and 
Ogata discussed a competition between the triplet p x - 
wave and a singlet (4,2 _ a 2 -wave, we find here that the 
real competitor (the most stable singlet state) is unex- 
pectedly an extended s-wave rather than d x 2_ y 2--w&ve. 
This is understood as follows. Since K(Q) > is re- 
pulsive, s (k)</> s (k + Q) < has to be satisfied. For the 
d x 2_ y 2-wa,ve, the k, k' = k + Q on the Fermi surface that 
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FIG. 2. Pair scatterings across the Fermi surface which 
favor (O) °r unfavor (x) p x -wa,ve pairing (a), those con- 
tributing to d x 2_ y i (b) or to extended a pairings (c). The 
white(grey) areas represent positive(negative) (f>. 
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FIG. 3. Eigenvalues of the Eliashberg equation for triplet 
and singlet pairings as a function of temperature for (7 = 5, 
t' = 0.3, a = 0.8. 




along the Fermi surface is displayed in Fig.|4|, which 
has a crib shape— with dips at the corners around k = 
(±27r/3, ±27r/3)B The dip arises because 4>t in each of 
the p x and p v channels has already sharp drops at the 
corner of the Fermi surface. Thus, in phase-insensitive 
experiments, the gap function obtained here may look 
like a two-dimensional /-wave pairing with nodes along 
k x = ik v , since a dip and a node are indistinguishable 
when the temperature is greater than the dip. Further 
study on this point is under way. 

Discussions with Yuji Matsuda are gratefully acknowl- 
edged. This work is in part funded by the Grant-in- Aid 
for Scientific Research from the Ministry of Education of 
Japan. 
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FIG. 4. |0(k)| (eq.(0)) for U = 5, t' = 0.3, a = 0.8, and 
T = 0.02. The inset shows how the two functions are com- 
bined. 

satisfy k - k' - (±2tt/3, ±2tt/3) (mod(27r, 2tt)) and 
a (k)0 a (k') < are only k,k' - (±2tt/3, 0), (0, ±2tt/3) 
(dashed arrows in Fig.||(b)). By contrast, all the k, k"s 
on the Fermi surface with k — k' ~ (±2tt/3, ±27r/3) con- 
tribute to the extended s-wave pairing (Fig.||(c)). 

The competition between the triplet and singlet is 
quantified by the eigenvalues of the Eliashbexg equation, 
shown in Fig.|| as functions of temperature.EJ It can be 
seen that for the value of a(= 0.8) adopted here, the 
triplet does dominate at low temperatures. The singlet 
extended-s, although having a greater magnitude of the 
pairing interaction, has the nodal line running in the 
vicinity of the Fermi surface (Fig.^(d)), which should be 
why this pairing is weaker. If we extrapolate A to T — > a 
finite transition temperature much smaller than O(0. Oli) 
is suggested. 

The above result and argument are for the p^-wave. 
Obviously, the p y -w&ve should be degenerate with it from 
symmetry. A complex linear combination of these two, 
Px J rip yi which breaks the time reversal symmetry, should 
be the true state below T c since the superconducting gap 
is maximized for that combination. 

The absolute value of the gap function for that linear 
combination on the Fermi surface, 

|0(k)| = [Mkx,ky) 2 + Mky,k x ) 2 ] 1/2 , (11) 
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The simultaneous enhancement of the triplet pairing and 
the spin fluctuations is reminiscent of the Tomonaga- 
Luttinger theory for purely ID systems, where one can 
show that the susceptibilities for the triplet pairing and 
the SDW may diverge simultaneously for T — > when the 



SU(2) symmetry in the interaction is broken. However, in 
the Tomonaga-Luttinger case, the triplet pairing dominates 
over the SDW only when the electron-electron interaction 
is attractive, in sharp contrast with the present case where 
the interaction is repulsive. The difference comes from a 
different origin of the broken SU(2) here, which is the 
anisotropy in the irreducible susceptibility (which arises at 
a single electron level) rather than the anisotropic electron- 
electron interaction. 

We note that, while the qualitative result does not depend 
crucially on t'(= 0.3 in this paper), the quantitative shape 
of the dip does. 



